Universality of electron correlations in conducting carbon nanotubes 



OO 
On 
On 

> 
o 

O 



i 

C 

o 
o 



> 



oo 

ON 
-i— > 

c3 



I 

C 

O 

o 



X 



Arkadi A. Odintsov 1,2 and Hideo Yoshioka 2 ' 3 
1 NEC Research Institute, 4 Independence Way, Princeton, New Jersey 08540 
2 Department of Applied Physics, Delft University of Technology, 2628 CJ Delft, The Netherlands. 
3 Department of Physics, Nagoya University, Nagoya 464-8602, Japan. 
(February 1, 2008) 

Effective low-energy Hamiltonian of interacting electrons in conducting single-wall carbon nan- 
otubes with arbitrary chirality is derived from the microscopic lattice model. The parameters of the 
Hamiltonian show very weak dependence on the chiral angle, which makes the low energy properties 
of conducting chiral nanotubes universal. The strongest Mott-like electron instability at half filling 
is investigated within the self-consistent harmonic approximation. The energy gaps occur in all 
modes of elementary excitations and estimate at 0.01 — 0.1 eV. 



PACS numbers: 71.10.Pm, 71.20.Tx, 72.80. Rj 

Single wall carbon nanotubes (SWNTs) are linear 
macromolecules whose individual-properties can be stud- 
ied by methods of nanophysicsEl. Recent demonstrgu 
tion of electron transport through singlea and multipletl 
SWNTs has been followed by remarkable observa- 
tions of atomic structurecm, one-dimensional van Hove 
singularities!], standing electron wavesa and, possibly, 
electron correlationsQ in these systems. Moreover, the 
first prototype of a functional device - the nanotube field 
effect transistor working at room temperature - has been 
fabricated recently£3. 

Structurally uniform SWNTs can be characterized by 
the wrapping vector w = Nidi + A^a^ given by the 
linear combination of primitive lattice vectors a± = 
(±l,v / 3)a/2, with a w 0.246 nm (Fig. 1). ft is natu- 
ral to separate non-chiral armchair (Ni = A2) and zig- 
zag (Ni = —A/g) nanotubes from their chiral counter- 
parts. Recent scanning tunneling microscopy studyo has 
revealed that individual SWNTs are generally chiral. Ac- 
cording to the single-particle model, the nanotubes_with 
N% — N2 = mod 3 have gapless energy spectrumQ and 
are therefore conducting; otherwise, the energy spectrum 
is gapped and SWNTs are insulating. Therefore, on the 
level of non-interacting electrons, physical properties of 
SWNTs are determined by their geometry. 

The Coulomb interaction in one-dimensional SWNTs 
should result in a variety of correlation effects due to 
the non-Fermi liquid ground state of the system. In 
particular, metallic armchair SWN3^|_ace. predicted to 
be Mott insulating at half-nllingpji3,y,y. Upon dop- 
ing the nanotubes become conducting but still display 
density wave instabilities in three modes of elementary 
excitations with neutral total chargeOjia. 

Experimental observation of electron correlations still 
remains challenging since their signatures are usually 
masked by charging effects. Two recent transport 
spectroscopy experiments on .an individual conducting 
SWNTd and a rope of SWNTsEiproduced contradicting 
results. The data by Tans et. alB assumes spin polarized 
tunneling into a nanotube, which in turn suggests the 
interpretation in terms of electron copulations. On the 
other hand, the data by Cobden et. al.Lj fits the constant 
interaction model remarkably well and shows no signa- 
tures of exotic correlation effects. Since the atomic struc- 
ture of particular SWNTs studied in these experiments 



is not known, it might be appealing to interpret the dis- 
crepancy in the results in terms of geometry-dependent 
many-particle properties of conducting SWNTs. 

Unfortunately, consistent theory of interacting elec- 
trons in chiral SWNTs is lacking, despite such a 
theory has., bp p .T recently developed for armchair 
nanotubesE£l,Ej,EJjl3. In this work we establish effective 
low-energy model for conducting chiral SWNTs and eval- 
uate its parameters (scattering amplitudes) from the mi- 
croscopic theory. We found very weak dependence of the 
dominant scattering amplitudes on the chiral angle. This 
allows us to introduce universal low-energy Hamiltonian 
of conducting SWNTs. According to the results of the 
renormalization group analysistj the strongest Mott-like 
electron instability occurs at half-filling. We investigate 
this instability using the self-consistent harmonic approx- 
imation. Substantial energy gaps are found in all modes 
of elementary excitations. The conditions for experimen- 
tal observation of the gaps are briefly discussed. 

We start from the standard kinetic term Hk of the 
tighti-binding Hamiltonian for p z electrons of a graphite 
sheeJla, 



Hk — 



J2{t(k)al jS (k)a + , s (k) + h.c.} 



(1) 



Here a PjS (fc) are the Fermi operators for electrons at the 
sublattice p = ± (Fig. 1) with the spin s = ± and 
the wavevector k = (k x ,k y ). The matrix elements are 
given by = -i( e - ifc v°/v^ + 2e ife « a / 2v/3 cos k x a/2), t 
being the hopping amplitude between neighboring atoms. 
The eigenvalues of the Hamiltonian vanish at two Fermi 
points of the Brillouin zone, aK with a = ± and K = 
(K,0), K = 47r/3a. 

We consider conducting chiral (ATi, A2) SWNT of the 
radius R = (a/2n)y / Nf + NiN 2 + AT| whose axis x' 
forms the angle x = arctan[(A^ 2 - A r i)/v / 3(A^i + N 2 )} 
with the direction of chains of carbon atoms (x axis 
in Fig. 1). Expanding Eq. (|l|) near the Fermi points 
to the lowest order in q = k — aK = q(cosx,sinx) 
and introducing slowly varying Fermi fields ip pas (x') = 
i_1/2 E F 2»/£ e. lqx 'a pa s{q + aK), we obtain, 
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FIG. 1. Graphite lattice consists of two atomic sublattices 
p — +, — denoted by filled and open circles. SWNT at the 
angle x to x axis can be formed by wrapping the graphite 
sheet along w = OO' vector. 



H k = -j^ae" 1 "^ / dx'^ pas d x ^. pas , (2) 

pas 

with the Fermi velocity v = V3ta/2 rs 8.1 X 10 5 m/s. The 
kinetic term can be diagonalized by the unitary transfor- 
mation 

^^^e-^^H 1 ?^ (3) 

v 2 r=± 

to the basis ip ra s of left- (r = — ) and right- (r = +) 
movers. 

The Coulomb interaction has the form 
H in t = \ V pA 2 <xK) 

pp' ,{ati} ,ss f 

x J dx'^l aiS 'tpl, a2S ,'tp p/a3Sl 'tp paiS , (4) 
with the matrix elements 
V pp '(2aK) = - U(r p - f p >) exp[-2iaK(r p - ry)], 

(5) 

corresponding to the amplitudes of intra- (p = p') and 
inter- (ja = — p') sublattice forward (a = 0) and back- 
ward (a = ±1) scatteringLJ (the sum is taken over 
the nodes of the sublattice p of SWNT). Here U(r) = 

e 2 /{n^Ja 2 + {x') 2 + 4i? 2 sm 2 (y'/2R)} is the Coulomb in- 
teraction with a short-distance cutoff ao ~ a, p = 
inR/^/3a 2 is a linear density of sublattice nodes along 
SWNT, and a = {a x - a 4 )/2 = (a 3 - a 2 )/2. We will 
choose the parameter ao from the requirement that the 
on-site interaction in the original tight-binding model 
corresponds to the difference between the ionization po- 
tential and electron affinity of sp 2 hybridized carbonE3. 
This procedure gives a = 0.526a. 



The forward scattering part Hp [terms with a — in 
Eq. of the Hamiltonian Hi nt can be separated into 
the Luttinger model-like term H p and the term Hf re- 
lated to the difference AV(0) — V pp (0)~ V p - P (0) between 
intra- and intersublattice amplitudes, 

H P = ^ / d X >p 2 ( X >), (6) 

H f = ~ AV 2 Q ^ J ^'^pas^-pcL'A-P^s'^pai-, (7) 

paa' ss' 

where p = J2 P as ^pas^pas is the total electron density. 
The backscattering Hamiltonian Hb [terms with a = ±1 
in Eq. ([|)] can be subdivided into the intrasublattice 

(p — p') and intersublattice (p — —p') parts. 

The dominant contribution to the forward scattering 
amplitudes V pp i (0) comes from the long range component 
of the Coulomb interaction, V pp (0) = (2e 2 / k) ln(R s /R), 
where R s — min(L, D) characterizes the screening of 
the interaction due to a finite length L of the SWNT 
and/or the presence of metallic electrodes at a dis- 
tance The forward scattering differential part 
A]/(0) and the intrasublattice backscattering amplitude 
V PP (2K) can be estimated from Eq. (||) as follows, 
AV(0), V PP (2K) ~ ae 2 /nR (for a ~ a). Despite the am- 
plitudes AV(0), V PP (2K) are much smaller than V pp (0), 
they produce essentially non-Luttinger terms in the low- 
energy Hamiltonian which will be important in further 
analysis. 

We evaluated the matrix elements (||) numerically for 
chiral SWNTs with radiuses R in the range 2R/a = 4 — 7 
(2R/a = 5.5 for (10,10) SWNTs). We found that dimen- 
sionless amplitudes 2ttkR[AV(0), V pp (2K)]/ae 2 show 
very weak dependence on the radius of SWNT and its 
chiral angle (see Table 1). The results are sensitive to 
the value of the cutoff parameter ao- 

The intersublattice backscattering am- 
plitude Vp- P (2K) is almost three orders of magnitude 
smaller than AV(0), V PP (2K). This is due to the C 3 
symmetry of a graphite lattice, which leads to an exact 
cancellation of the terms (||) contributing to V P - P (2K) in 
the case of a plane graphite sheet. The matrix elements 
V P - P (2K ) are generally complex due to asymmetry of ef- 
fective ID intersublattice interaction potential (the ma- 
trix elements are real for symmetric zig-zag and armchair 
SWNTs). Let us note that after the unitary transforma- 
tion (||) of the Hamiltonian H = Hk + Hi nt , the chiral 
angle x enters mly to the intersublattice backscattering 
matrix elementstll. Due to the smallness of these matrix 
elements, the low-energy properties of chiral SWNTs are 
expected to be virtually independent of the chiral angle. 



ao/a 


AV(0) 


V PP (2K) 


\V P - P (2K)\ 


0.4 

0.526 

0.7 


0.44265 - 0.44274 
0.17378 - 0.17395 
0.04880 - 0.04895 


0.97060 - 0.97095 
0.53549 - 0.53561 
0.24778 - 0.24797 


0.6 - 2.2 x 10" a 
0.5 - 1.6 x 10" 3 
0.3-1.5 x 10" 3 



TABLE I. Scattering amplitudes AV(0), V PP (2K), 
Vp- P (2K) in units ae 2 /2-kkR for all chiral SWNTs with ra- 
diuses 7? in the range 2R/a = 4 — 7. 
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Neglecting the intersublattice backscattering we arrive 
to the universal low-energy model of conducting SWNTs 
given by the Hamiltonian H = H k + H F + H^— The 



latter can be bosonized along the lines of Refs.EjO. We 
introduce bosonic representation of the Fermi fields, 



^Pras 



Vras 



: exp 



vrq F x' + — {6 as + r(j) as } 



(8) 



and decompose the phase variables 9 as , <f> aa into symmet- 
ric 5 = + and antisymmetric S = — modes of the charge p 
and spin a excitations, 6 as = 9 p+ + s9 a+ + a9 p - + as9 a - 
and 4>as — 4> P ++s(j) cr+ +a(j) p -+as4> cr - . The bosonic fields 
satisfy the commutation relation, [6j${xi), <f>j>s'(^2)] — 
i(7r/2)sign(:ri — xaWjj' $55' ■ The Majorana fermions 
i] ras are introduced^ to ensure correct anticommutation 
rules for different species r, a, s of electrons, and satisfy 
[Vras,r)r>a's']+ = 26 rr >5 aa >6 ss > . The quantity q F = im/A 
is related to the deviation n of the average electron den- 
sity from half-filling, and a ~ a is the standard ultraviolet 
cutoff. 

The universal low-energy Hamiltonian of conducting 
SWNTs has the following bosonized form, 



# = E E? Ux'{Kj 5 l {d x ,e jS f + K j5 {d x 4iS?} 

A — „ X L •* 

dx'{[AV{0)-V pp {2K)] 



j—p,cr S=± 
1 

[cos(4g F ir' + 29 p+ ) cos 29 a+ — cos 29 p _ cos 29^-} 
-AV(0) cos(4q F .x' + 29 p+ ) cos 26> p _ 
+AV{0) cos(4g F .x' + 29 p+ ) cos 20 CT _ 
-AV(0) cos2(9 CT+ cos20 p _ 
+AV{0) cos26> CT+ cos20 CT _ 
-Vpp{2K) cos(4g F x' + 29 p+ ) cos 2<j) (T - 
+V PP (2K) cos 29 a+ cos 20 CT _ 
+V PP (2K ) cos 26 p - cos 20 CT _ 

+1^(21?) cos 20a- cos 2^-}, (9) 



f J( 5 = v<J AjgBjs and i?^ = ^Bjg/Ajs being the veloci- 
ties of excitations and exponents for the modes j, 6. The 
parameters Ajg, Bjs are given by 



A„ 



A„ 



4^(0) AV{Q) V PP (2K) 

ttv Attv 2ttv 

AV(0) S V PP {2K) 



B v s — 1 



Am 
AV(0) 

47TW 



2ttv 



(10) 

(11) 

(12) 



with V-(O) = [^p(O) + V^_ p (0)]/2. 

Since V(0)/w ~ ln(i? s /i?) > 1 and AV(0)/v, 
^,p(2i4T)/w ~ a/i? <C 1, the renormalization of the pa- 
rameters Kjg, vjs by the Coulomb interaction is the 
strongest in p+ mode. Assuming k — 1.4E3 and R s = 100 
nm we obtain K p+ ~ 0.2.-u p+ ~ v/K p+ for generic 
SWNTs with i? - 0.7 nm£j. The interaction in the 
other modes is weak: Vjs = V, Kjg = 1, up to a fac- 
tor 1 + 0(a/R). 



The renormalization group analysis of armchair 
SWNTs with longjABge Coulomb interaction has been 
performed in Refs£jji3. The modification of the param- 
eters of the Hamiltonian (g) by the neglected small tepa 
Vp- P (2K) should not change the results qualitativelyea. 
The most relevant perturbation is the umklapp scat- 
tering at half-filling. In this case the non-linear terms 
of the Hamiltonian (Q) , which do not contain cos 29 a - 
scale to strong coupling and the phases 9 p+ , 9 a+ , 6> p _ , 

get locked at (9 p +\ 9^ , 9 p m \ ^ ) = (0,0,0,0) or 
(tt/2, 7r/2, 7t/2, 7r/2). Therefore, the ground state of half- 
filled SWNT is the Mott insulator with all kinds of the 
excitations gapped. 

To estimate the gaps quantitatively, we will employ 
the self-consistent harmonic approximation which follows 
from Feynman's variational principleE3. We consider trial 
harmonic Hamiltonian of the form: 



H ° = E If f dx ' i K l5 l \{d xl e jS f + (1 - & ja 8 S -)<$ s P jS ] 

+ K ]S [(d xl <p jS ) 2 + 6 ja 6 S -q] s cf> 2 jS }}, (13) 

qjS being variational parameters. By minimizing the up=. 
per estimate for the Free energy F* = F a + (H - Ho)(fB 
one obtains the following self-consistent equations, 



2K n 



■-(>+ 



-c p+ {[V pp {2K)-AV(0)}c, 



+AV(0)c p ^+V pp (2K)d a _}, 



2K n 



%- = 



-c p _ {Ay(0)c p+ + AV(0)c, 



V pp {2K)d a .} . 



2K, 



Trd 2 v eT + 



{[V pp (2K)-AV(0)]c p+ 



+Al/(0)c p _ - V pp (2K)d a -} , 
2 



(14) 
(15) 
(16) 



Qa- = 



na 2 K„-V„- 



d a _V pp {2K) {c p+ - c a+ - c p _} , (17) 



where c 



(cos 29 



29 { ™\ 1 hq 3 s) K ^ , d a _ 



(cos20 cr _) o = cos2</)^ ("faqa-) 1 ^"- , (...) denotes av- 
eraging with respect to the trial Hamiltonian (llH), and 
7 ~ 0.890 for the exponential ultraviolet cutoff. Note 
that (cos 2# (T _) = 0, so that only the terms of the Hamil- 
tonian (|j) which scale to the strong coupling contribute 
toEqs. !!)-©■ 

In the limiting case of interest, |Ay(0)|, |y pp (2i : C)| -C v 
and K p+ 1, the solution of Eqs. ([l^)-([L7|) can be found 
in a closed form, giving rise to the following estimates for 
the gaps Ajs — Vjsqjs in the energy spectra, 



A P + = 
A P - = 
A CT+ = 
A„_ = 



2l 2 V f 



1/(1-A' P+ ) 



A 



ja 

|AV(0)| 
V P + 

\V PP (2K )- AV(0)\ 

V P+ 
\V PP (2K )\ 



A 



Vc 



V+ 



(18) 
(19) 
(20) 
(21) 
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FIG. 2. The energy gaps Ajg for the modes p+, a—, a+, 
p— at ao = 0.526a (lines marked by crosses, from top to 
bottom) and for the mode p+ at ao = 0.4a (triangles) and at 
do — 0.7a (squares). The energy is in units Tiv/a ~ 2.16 eV 
for a = a. 

with V p+ = {[AV(0)-V PP (2K)] 2 + [AV(0)] 2 

1 /2 

+ [V PP (2K)) 2 } 1 . In the above expressions we used the 
approximation, vq/v p + = K p + and vo/vjs — Kjs = 1 
for <7± and p— modes. The formulae (p^)-(pl]) indi- 
cate that the largest gap occurs in p+ mode, albeit all 
four gaps are of the same order for realistic values of 
the matrix elements (see Table 1). The gaps decrease 
as A j5 (x (l/R) 1 ^ 1 -*^ ~ (1/i?) 5 / 4 with the tube ra- 
dius. This should be contrasted to the 1/i? dependence 
of wide semiconductor gaps and 1 / R 2 dependence of nar- 
row deformation induced gapsa expected from the single 
particle picture. 

In Figure 2 we present numerical results for the gaps 
Ajs for the cutoff parameter ao = 0.526a. The data for 
somewhat larger and somewhat smaller values of ao indi- 
cate possible variation of the gap A p+ due to uncertainty 
in the short distance cutoff of the Coulomb interaction. 
The gaps can be loosely estimated at Aj$ ~ 0.01 — 0.1 eV 
for typical SWNTs with R ~ 0.7 nm. Due to the gaps in 
the spectrum of bosonic elementary excitations, the elec- 
tronic density of states should disappear in the subgap 
region and display features at the gap frequencies and 
their harmonics. Both signatures should be observable 
by means of the tunneling spectroscopy. 

Why hayep- .the gaps not been observed in the 
experimentaJ'ESu? This might be due to the effect of 
metallic electrodes. The difference in the workfunctions 
of the electrodes (Au, Pt) and the nanotube results in 
a downward shift of the Fermi level of the nanotube by 
a few tenths of an eVa. This causes substantial devia- 
tion An = Aqp/n ~ 1 nm" 1 of the electron density in 
SWNT from half filling, at least in the vicinity of the 
electrodes. Therefore, we expect the gap features to be 
observable in the layouts with well separated (to a dis- 
tance d 3> TlvfIA p+ ) source and drain contacts. The 
piece of nanotube between them should be well isolated 
from any conductor. 



In conclusion, we have developed effective low-energy 
theory of conducting chiral SWNTs with the long-range 
Coulomb interaction. The many-particle properties of 
SWNTs are found to be virtually independent of the chi- 
ral angle. The universal Hamiltonian (||) of conducting 
SWNTs is introduced. The Mott-like energy gaps in the 
range of 0.01 — 0.1 eV should be observable at half filling. 

The authors would like to thank B.L. Altshuler, 
G.E.W. Bauer, R. Egger, Yu.V. Nazarov, and N. 
Wingreen for stimulating discussions. The financial sup- 
port of the Royal Dutch Academy of Sciences (KNAW) is 
gratefully acknowledged. One of us (A.O.) acknowledges 
the kind hospitality at the NEC Research Institute. 



1 A. Thess et. al, Science, 273, 483 (1996). 

2 S.J. Tans, et. al, Nature 386, 474, (1997). 

3 M. Bockrath, et. al, Science, 275, 1922 (1997). 

4 J.W.G. Wildoer, et. al, Nature 391, 59 (1998). 

5 We do not discuss here the minieaps due to the curvature 
of the nanotube surface, see RefJj. 

6 C.L. Kane, E.J. Mele, Phys. Rev. Lett. 78, 1932 (1997). 

7 W. Clauss, D. J. Bergeron, A. T. Johnson, Phys. Rev. B 
58, R4266 (1998). 

8 L.C. Venema, et. al, submitted to Science. 

9 S.J. Tans, M.H. Devoret, R.J.A. Groeneveld, and C. 
Dekker, Nature 394, 761 (1998). 

10 S.J. Tans, A.R.M. Verschueren, and C. Dekker, Nature 
393, 49 (1998). 

11 L. Balents and M. P. A. Fisher, Phys. Rev. B 55, R11973 
(1997). 

12 Yu. A. Krotov, D.-H. Lee, and Steven G. Louie, Phys. Rev. 
Lett. 78, 4245 (1997). 

13 C. Kane, L. Balents and M. P. A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 

14 H. Yoshioka and A. A. Odintsov, to be published in Phys. 



Rev. Lett., preprint 
R. Egger and A. O 
(1997). 

D.H. Cobden, et. al, Ph; 
Drawbacks of the analysi; 
P. R. Wallace, Phys. Rev 



ond-mat/9805106 



Gogolin, Phys. Rev. Lett. 79, 
Rev. Lett. 81, 681 (A 



5082 



are discussed in Ref£fl. 
71, 622 (1947). 

19 This definition is different from conventional one since it 
involves the Fermi points rather than the directions of mo- 
tion of scatterred electrons. 

' see e.g. E. Moore, B. Gherman, and D. Yaron, J. Chem. 
Phys. 106, 4216 (1997). 

The transformed matrix elements V p - p (2aK)e 2lapx are in- 
variant under 7rn/3 rotations in chiral angle \ n — \o + 
nn/3, which correspond to equivalent representations of a 
nanotube. 

; We have checked this explicitly for (10,10) SWNTs at half 
filling and away from it. 

' R.P. Feynman, Statistical mechanics, Reading, Mass., 
W.A. Benjamin, Inc., 1972. 



4 



